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Abstract 

We use an off-lattice discretization of fractional Brownian motion 
and a Metropolis Algorithm to determine the asymptotic scaling of 
this discretized fractional Brownian motion under the influence of an 
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excluded volume as in the Edwards and Domb-Joyce models. We find 
a good agreement between the Flory index describing the scaling of 
end-to-end length with a mean field formula proposed earlier for this 
class of models. 


1 Introduction 


1.1 Background 

Random paths and in particular (weakly) self-avoiding paths have 
been studied intensively in physics as well as in probability theory. 
In physics they play a crucial role in the modeling of chain polymers 
with ’’excluded volume effect”. In these models self-intersections are 
penalized; in the continuum case this is the Edwards model I IEdw65j . 
in the discrete case one speaks of the Domb-Joyce model [ DJ72l |. The 
penalization factor in the Edwards model is introduced through a 
modification of the Wiener measure go by a Gibbs factor 


d/x 




dt5{B{r) — B(t)) 


dg 0 ) 


where N, g > 0, 5 is the d-dimensional Dirac delta function, B denotes 
a version of a d-dimensional Brownian motion and 


rN 


Z = El exp —g / dr dt5(B(r) — B(t)) 


We set 

rN rN 

L= / dr dt 5 (B(t) — B(t)), 

Jo Jo 

which is known as the ’’self-intersection local time” of a Brownian 
motion. 

The one-dimensional weakly self-avoiding random walk is well un¬ 
derstood and proofs of mathematical rigor concerning the scaling ex¬ 
ponent can be found in |GdH93) . 

For fractional weakly self-avoiding paths the contrary is the case. 
This is due to the lack of Markov and martingale properties, which 
excludes many techniques of stochastic analysis. 
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1.2 Fractional Model 


Fractional Brownian motion (fBm) has been suggested as a more gen¬ 
eral model for chain polymers, see e.g. jHMR.1 3] . 

FBm on M d , d > 1, with ’’Hurst parameter” H e (0,1) is a d- 
tuple of independent centered Gaussian processes iBHOZOfl IMis08| 
B h = {Bf 1 : t> 0} with covariance function 


E(#f B?) = * ( t 2H + t 2H - 1 1- t\ 2H ) . (1.1) 


One sees that, for H = 1/2, one recovers ordinary d-dimensional 
Brownian motion. From (1.1) one further verifies stationary incre¬ 
ments: 


E 


[(B H (t)-B H (r)) 2 ) =\t-r\ 


2 H 


and since the increments are centered Gaussians, we have more gen¬ 
erally 


E (( B H (t) - B H {r)) 2k ^j = c k 1 1 - r\ 2Hk , c k > 0, 

which, by the Kolmogorov-Chentsov theorem implies the continuity 
of sample paths. Increments are uncorrelated for H = 1/2, corre¬ 
lated (’’persistent” paths, or stiffer polymers) for H > 1/2, and anti¬ 
correlated (curlier than Brownian motion) for small Hurst indices 
H < 1/2. Because of these properties fractional Brownian motion 
paths have been proposed as generalized chain polymer models by 
|BC95| A rigorous mathematical extension of the Edwards model was 
recently shown jGOdSSlTj to exist for a limited range of Hurst param¬ 
eters H and dimensions d. This limitation is due to the singular nature 
of the self-intersection local time and will not be shared by discretiza¬ 
tions which are thought to be in the same universality class of asymp¬ 
totic scalings. The characteristic observable of this general class of 
chain polymer models is the average mean-squared end-to-end length 
R 2 and its asymptotic scaling behavior as the number of monomers 
increases. In this paper we intend a first exploration of the scaling 
behavior for fBm-based models. We should emphasize that we do not 
yet aim for high precision results in this first study. In this exploratory 
phase we have restricted ourselves to the one-dimensional case. There 
the mean end-to-end length of paths grows linearly for self-avoiding 
random walks |vdHdHK97) . as is to be expected intuitively. For frac¬ 
tional paths with H < \ on the other hand the situation is far from 
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being obvious: the repulsive excluded-volume dynamics is balanced 
by long-range attractive force acting along fractional Brownian trails 
for small Hurst index. These models might be considered unphysical. 
However, as in the Brownian case, there exists a recursion formula 
}KF78| . which would predict the Flory index in higher dimensional 
cases from the one-dimensional one, mm- 


2 The Algorithm 

2.1 Discretization and Metropolis Algorithm 

There are various methods to simulate fBm paths. Exact meth¬ 
ods such as the ones of Hosking |Hos84| . the Cholesky method, see 
e.g. [ AG07 ] or the method of Wood and Chan }CW97i |. which uses roots 
of the covariance to obtain a fractional path from a standard normal 
distribution are not usable here, since it is not obvious how to include 
directly the excluded volume. The same holds for approximative meth¬ 
ods using spectral densities, the Paxson method {Pax97j . which uses 
Fourier transform or the Decreusefond-Lavaud method |DL96| which 
is using the kernel representation of fBm. 

A more natural way to implement a discrete fractional walk with self- 
repellence is to use off-lattice discretization and Monte Carlo meth¬ 
ods based on a Metropolis algorithm as in |BGZ00] for the Brownian 
case and in lll.MR 1.1 for the fractional Brownian motion. In the our 
algorithm the Metropolis routine is used to update ’’bonds”, i.e. in¬ 
crements of the path. Given a Hurst parameter H E (0,1), we define 
the positions of a walk with N points by 

Xj = B H {j), j = 0, 1, 

and N — 1 ’’bond vectors” 

y j = B H (j + l)-B H (j). 

These are centered Gaussian, so one obtains their probability density 
by inverting the covariance Aji = E (ijj yi ). 

6o(y) = , 1 = e x P ( ) = C ex P \-\{y,U Q y) 

yEWETh v 2 ) v 2 
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To obtain a one-dimensional weakly self-avoiding fractional random 
walk we have to penalize self-crossings. To discretize the self-intersection 
local time: 

pN rN r 

L= dr dt5 [B h (t) — B H (t)) = / duL 2 u , 

Jo Jo Jr 


with 


rN 


L u = dt 5 (B h ( t) — u) , u G M, 

Jo 

we decompose M into intervals I n of equal length l and replace L u 
the number of positions Xk that fall into the interval with number 


by 

n: 


L n — \^j £ In}■ 

Likewise we replace the self-intersection local time by 

^ = E L «> 

n 

hence the unnormalized probability density of conformations becomes 

p{x) ~ exp (-^(y,H 0 y) - g(L(y)) 

The well-known Metropolis algorithm, see e.g. |AG07j in our case with 
the energy 

E = [y.'Hoy) - gL{y), ( 2 . 1 ) 

does not require normalization of the probability density. In the above 
formula y G is the vector of increments and g > 0 is the coupling 

constant of the excluded-volume term L. 



2.2 Qualitative Description of the Algorithm 

For a fixed polymer length N we start from a random initial configura¬ 
tion, randomly chosen bond vectors are updated and subjected to the 
standard Metropolis routine. After an initial relaxation phase Nr, the 
end-to-end lengths are sampled, not after individual updates but in 
intervals large enough to suppress correlations. Finally the average of 
the sample is determined. The process then is repeated for increasing 
N to determine the scaling of the end-to-end length with growing N. 
The error analysis was based on standard analysis of 30 independent 
runs with in particular individual random seeds. 
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2.3 Parameter Choice 


For the number of monomers N our choice of the range between 200 
and 600 was dictated by two considerations. N should be large enough 
to approximate the asymptotic regime. The upper limit was chosen 
such that equilibrium could be reached within a few days of computing 
time. The interval length / in the definition of incidences above was 
chosen to be l = 1. 



Coupling constant g 

Figure 1: Dependence of the scaling exponent from the penalizing strength 
g (for r — 3 and H = 0.5, N = 100,120,160 and 200, s = 25 x 10 6 updates 
and Nr = 5 x 10 6 relaxation time with 30 repetitions). 

Recall that g is the coupling constant or strength of excluded vol¬ 
ume effect. For g = 0 the walks have no penalization and will scale as 
a fBm-path, i.e. the scaling exponent will be v = H. For small posi¬ 
tive g the scaling index u rises sharply, see Figure 1 for the Brownian 
case. It overshoots the theoretical value and then seemingly drops 
below it for larger g. Scrutiny of the end-to-end length in this regime 
reveals the cause of this. As g increases there appears an increasing 
number of ’’outliers” with much shorter end-to-end length, i.e. confor¬ 
mations that did not unfold to equilibrium conformations during the 
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relaxation and sampling period. To amend this problem we did not 
perform relaxations with a constant g, instead we slowly increased to 
coupling constant from zero, with the result that the outliers disap¬ 
peared and the scaling index stabilized. See the black line in Figure 
1. Monte Carlo updates of bond vectors are chosen from a sampling 



Sampling Window r 


Figure 2: End-to-end length as a function of sampling window r (for g = 
2 ,H = 0.5, N = 200,250,300, s = 25 x 10 6 and Nr = 5 x 10 6 with 30 
repetitions). 

window [—r, r]. As a function of r the computed end-to-end length 
will saturate within the sampling period when r is large enough, see 
Figure 2. This led us to choose r = 4. Larger sampling windows would 
increase the rate of rejects in the Metropolis algorithm. 

An analogous analysis for Hurst parameters H < \ produces the 
same results for suitable g and r. 

One observes that the relaxation of end-to end length is some or- 


7 





ders of magnitudes slower than that of the two energies in equation 
(2.1), [Cabl4j . After some studies we chose 5 x 10 6 relaxation updates 
Nr as a matter of computational convenience; which proved to be suf¬ 
ficient when we modified the relaxation protocol as described above 
to facilitate unfolding. The number of updates in the sampling phase 
was s = 25 x 10 6 , and sampling was performed every 10000 updates. 


3 Results and Conclusions 


3.1 The Flory index 


While there are not yet any mathematical proofs it is generally ex¬ 
pected that for large N the end-to-end length of polymers will scale ac¬ 
cording to R 2 ~ N 2u when they are modeled by (weakly) self-avoiding 
random paths. Mean field arguments - notoriously unsound and yet 
successful - would suggest the following formula for the scaling index 
in the fractional case 


v = 


2H + 2 
d + 2 ’ 


(3.1) 


with N 1 generalizing the famous Flory formula [ Flo 53] for H=\. 
In Figure 3 the prediction in (3.1) is compared to our numerical find¬ 
ings. 
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